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In the frame work of a nuclear overlap model, we estimate the number of nucleon and quark partic- 
ipants in nucleus-nucleus collisions. At RHIC and LHC energies, the number of nucleons-normalized 
charged particle and transverse energy density in pseudorapidity, which shows a monotonic rise with 
centrality, turns out to be an almost centrality independent scaling behavior when normalized to the 
number of participant quarks. This indicates partonic activities both in RHIC and LHC heavy-ion 
collisions. We study the charged particle and transverse energy production mechanism from AGS, 
SPS, RHIC to LHC energies in the framework of nucleon and quark participants as well. A universal 
function which is a combination of logarithmic and power-law, describes the charged particle and 
transverse energy production both at nucleon and quark participant level quite well for the whole 
range of collision energy under discussion. A detailed energy dependent production mechanisms 
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are discussed both for nucleonic and partonic level. We give the predictions for 

barometric observable, -^J- for Pb+Pb collisions at Jsnn = 5.5 TeV. A universal law of particle 

production from elementary to heavy-ion collisions is proposed. 

PACS numbers: 25.75.Ag, 25.75.Dw 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD), the theory of 
strong interaction, predicts a phase transition from nor- 
mal hadronic matter to a deconfined state of matter 
called Quark Gluon Plasma (QGP) [IH3. Such a phase 
transition occurs at very high temperature and/or energy 
density. It is possible to create these extreme conditions 
of high temperature and/or energy density in the labora- 
tory by colliding heavy-ions at ultra-rclativistic speeds. 
The properties of this new phase of partonic matter is 
yet to be fully understood. The transverse energy and 
charged particle multiplicity are two key observables in 
characterizing the bulk properties of the matter created 
in heavy-ion collisions. The dependence of charged mul- 
tiplicity and transverse energy on the collision geometry 
i.e. on centrality and collision energy is of paramount 
importance to understand the particle production mech- 
anism. This may give an insight to the partonic phase 
that might be created in heavy-ion collisions. 

In this paper, we calculate the number of nucleon 
and quark participants in a nuclear overlap model in a 
way done previously to study the enhancement of multi- 
strange baryons in heavy-ion collisions [J. We study 
the centrality and collision energy dependence of charged 
particle and transverse energy production. After finding 
out a universal law of particle production, we give the 
prediction for —^, ^§jf- and the barometric observable, 

in = 5.5 



di] 



SS s ® for Pb+Pb collisions at 

TeV. 
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II. CALCULATION OF THE NUMBER OF 
PARTICIPANTS 



Here, both nuclei and nucleons are considered as su- 
perposition of constituent or "dressed" quarks (partons 
or valons). There are three constituent quarks per nu- 
cleon. Baryons are composed of three and mesons are 
of two such quarks. The concept of constituent quarks 
is very well known [5J, |6[ and established in the realm 
of the discovery of constituent quark scaling of identi- 
fied particle elliptic flow at RHIC pj- The constituent 
quark approach is successful in explaining many fea- 
tures of hadron-hadron, hadron-nucleus |8| and nucleus- 
nucleus collisions. These include global properties like 
the charged particle and transverse energy density per 
participant pair [{J [lfj. QCD calculations support the 
presence of three objects of size 0.1-0.3 fm inside a nu- 
cleon [1JJ]. Further more it has been seen that nucleus- 
nucleus collisions and p+p collisions have similar initial 
states if the results are scaled by the number of con- 
stituent quark participants [12J |. These observations also 
indicate that the particle production is essentially con- 
trolled by number of constituent quarks pairs participat- 
ing in the collision. In a constituent quark picture, NN 
collision looks like a collision of two light nuclei with es- 
sentially one qq pair interacting in the collision, leaving 
other quarks as spectators. These quark spectators form 
hadrons in the nucleon fragmentation region with a part 
of the entire nucleon energy being used for the particle 
production (^/s qq ~ ^/snn/3). Contrary, in AA colli- 
sions, due to large nucleus size and the higher probability 
of quark-quark interaction from same projectile and dif- 
ferent target nucleons, more than one quark per nucleon 
take part in the interaction. 

The calculations of the mean number of nuclcon/quark 



participants are done in the following way. In the nu- 
clear overlap model, the mean number of participants, 
i.e. N^-part, in the collisions of a nucleus A and a nu- 
cleus B with impact parameter b is given by 
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where T(b) = J_ dz riA(\/b 2 + z 2 ) is the thickness func- 
tion, defined as the probability of having a nucleon- 
nucleon (NN) collision within the transverse area ele- 
ment db. [1 — aNNTA(b)/A] is the probability for a nu- 
cleon to pass through the nucleus without any collision. 
A and B are the mass numbers of two nuclei participat- 
ing in the collision process. We use the Woods-Saxon 
nuclear density profile, 



n A (r) = 



»o 



1 + exp[(r - R)/d] ' 



(2) 



with parameters, the normal nuclear density n = 
0.17 /to" 3 , the nuclear radius R = (1.12A 1 / 3 - 
0.86" 1 / 3 ) fm and the skin depth d = 0.54 fm. The 
inelastic nucleon-nucleon cross section, i.e. ctnn-i arc 
42 mb at ^/snn = 200 GeV and 64 mb at ^/snn = 
2.76 TeV [13| . In order to calculate the number of 
quark participants, iVq_ par t, in nucleus-nucleus collisions, 
the density for quarks inside the nucleus is changed to 
three times that of the nuclcon density (rig = 3n = 
0.51 fm~ 3 ). Instead of nucleon-nucleon cross section, 
quark-quark cross section is used which is 4.67 mb and 
7.1 mb at Js^ = 200 GeV and 2.76 TeV respectively 
0, i, OH El- For making predictions for Pb+Pb colli- 
sions at y^SNN = 5.5 TeV, we have used <tnn — 72 mb 
[15J and 8 mb for 7vN_ par t and iVq_ par t calculations, re- 
spectively. 



III. RESULTS AND DISCUSSION 

At lower center of mass energies, it has been found 
that the particle production scales with the number 
of participating nucleons, contrary to the case of high 
energies where hard processes dominate. Hard processes 
have much smaller cross-section than the soft processes. 
However, the number of binary collisions increase with 
increase in collision centrality faster than the number 
of participants, as a result the particle production per 
participant nucleon increases with increase in centrality. 
By using constituent quark approach, we are going to 
show how the particle production at higher energies 
depend on the participating quarks. For these studies, 
we have taken the centrality dependence of pseudo- 
rapidity density of charge particles and transverse energy 
from top RHIC energy to LHC V^nn = 2.76 TeV. To 
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FIG. 1: (Color online) Mean number of nucleon participants 
as a function of centrality in the overlap model (filled square) 
and from ALICE measurements (filled circles). The ratio of 
ALICE estimate and overlap model is shown in the lower 
panel, which shows a reasonably good agreement of both the 
estimations. 
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FIG. 2: (Color online) The ratio of iV q _ pa rt/AN-part as a func- 
tion of centrality for Pb+Pb collisions at ^/snn = 2.76 TeV 
from the overlap model calculations. This increases monoton- 
ically with collision centrality. 



study the constituent quarks dependence of charged 
particle and transverse energy production, we need 
to estimate the number of participant quarks, which 
has been done in the framework of nuclear overlap model. 

Since the quark participants are calculated in the 
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FIG. 3: (Color online) (a) 0.5 ./Visr_p art -normalized and (b) 
0.5 iV q _part-normalized d ch for Pb+Pb collisions at ^/snn = 
2.76 TeV for LHC and for Au+Au collisions at v /s^ = 200 
GeV for RHIC as a function of collision centrality. A scal- 
ing factor of 2.08 and 2.04 is multiplied with A^N-part and 
Afq-part-normalizcd RHIC data, respectively to make a direct 
comparison with the LHC data. Shown are the comparisons 
of IP-saturation model [13] and EKRT model [18| predictions 
for both the data sets. 



framework of nuclear overlap model, it is very essential 
to check how good is our estimation of number of par- 
ticipating nucleons in the collision. In order to do that, 
the mean number of participating nucleons, calculated in 
overlap model, are compared with the number estimated 
by the ALICE experiment [16[ . In Figure [TJ the mean 
number of participating nucleons are shown as a function 
of collision centrality. Solid circles represent the ALICE 
values and the solid squares represent the overlap model 
calculations. The lower panel of Figure [TJ represents the 
ratio of ALICE values and the overlap model values. It is 
clear from the figure that the overlap model calculations 
are in very good agreement with the ALICE calculations. 
We have then estimated the number of quark participants 
within the prescription described in the previous sections. 
The ratio of quark partipants and nucleon participants is 
shown in Figure [5] as a function of nucleon participants, 
which is a measure of collision centrality The ratio shows 



a sharp increase for iV 



N —part 



< 100 and then shows a 



type of linear monotonic rise going from peripheral to 
central collisions. This behavior is purely geometrical in 
nature. 

Figure [3] (a) describes the 0.5 ATN-part-normalizcd 
pscudorapidity density of charge particles for Pb+Pb col- 
lisions at y / snn = 2.76 TeV for LHC and for Au+Au 
collisions at -^/snn = 200 GeV for RHIC, as a function of 
collision centrality. A scaling factor of 2.08 has been mul- 
tiplied with the RHIC data to make a direct comparison 
with that of LHC data. Going from peripheral to central 
collisions the production of charged particles show a rise, 



when normalized with number of participant pairs. This 
rise in both the data sets are fitted with a function of the 
type: 
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where, A and n are fitting parameters. The value of 
n, is 0.18 ± 0.02 and 0.19 ± 0.02, for RHIC and LHC 
data sets, respectively. The value of A, is 2.72 ± 0.25 
and 2.63 ± 0.24 for RHIC and LHC data sets respec- 
tively. The expectations from IP-saturation model [131 
and the EKRT model [l8[ , which is based on initial state 
gluon saturation, are also shown for a direct comparison 
with data. In EKRT model the centrality dependence of 
charged particle multiplicity density in pseudorapidity is 
described by the function: 



^i h= o = c^i.i6(^r 2 (^)°- 40 

ar] 3 2 



(4) 



Figure 0(b) shows the 0.5 7Vq_ par t-normalized ^- £i - for 
Pb+Pb collisions at ^iiS" =2.76 TeV for LHC and for 
Au+Au collisions at ^/ifji = 200 GeV for RHIC as a 
function of collision centrality. A scaling factor of 2.04 
has been multiplied with the RHIC data set to make a 
direct comparison with the corresponding LHC data. A 
similar function as is given by Eqn. [3] has been fitted 
to both the data sets. The fitting parameter n, has the 
values -0.09 + 0.01 and -0.12 + 0.01 for RHIC and LHC 
data sets, respectively. The parameter A, has the val- 
ues 5.38 ± 0.40 and 6.2 ± 0.5 for RHIC and LHC data, 
respectively. With top peripheral collisions as excep- 



tions, the 0.5 N„ 



-normalized 



dJV c , 



both for RHIC 



'q-part "»"■ """""^ ^ 

and LHC data, show a centrality independent number of 
constituent quarks scaling behavior. The comparison of 
both the above figures reveals that the particle produc- 
tion is better described in terms of constituent quarks 
rather than participating nucleons. Hence, it makes sense 
to consider partons as the sources of particle production 
centers than nucleons. The nucleus-nucleus collisions at 
these energies are partonic in nature. 

The total transverse energy, which is measured on 
event- by- event basis, carries significant information re- 
garding the explosiveness of the reaction and on char- 
acterizing the global properties of the produced sys- 
tem. The estimation of transverse energy provides model 
based calculation of the initial energy density produced 
in the collision. This in turn gives information regard- 
ing a possible formation of a partonic phase, when com- 
pared with that of lattice QCD prediction of critical 
energy density for the deconfincment transition [3j. In 
nucleus-nucleus collisions, transverse energy is generated 
by the initial scattering of the partonic constituents of 
the incoming nuclei and by the rescattering of the pro- 
duced partons and hadrons. The transverse energy de- 
pends on the initial state of the collision and the viscos- 
ity of the partonic matter as it comes to the final state 
through interactions and expansion process. In Figure 0] 
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FIG. 4: (Color online) (a) 0.5 ./VN_p art -normalized and (b) 
0.5 iV q _p ar t-normalized ^^ for Pb+Pb collisions at ^/snn = 
2.76 TeV for LHC and for Au+Au collisions at v /s^ = 200 
GeV for RHIC as a function of collision centrality. A scaling 
factor of 2.8 is muplitiplied with the RHIC data in both the 
cases to make a direct comparison with LHC data. 



(a), 0.5 A r N_ pal . t -normalized ^^ for Pb+Pb collisions at 
^/snn = 2.76 TeV for LHC and for Au+Au collisions at 
-y/sNN = 200 GeV for RHIC are shown as a function of 
collision centrality. A scaling factor of 2.8 has been mul- 
tiplied with both the A^-part an d iVq_ par t-normalized 
RHIC data to make a direct comparison with the corre- 
sponding LHC data. Both the data sets show a functional 
trend which is governed by an equation, 
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N —part 
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= BN 



N—parf 



(5) 



where, B and n are fitting parameters. Here B has the 
dimension of energy. The values of B are, 2.9 ± 0.2 GeV 
and 5.8 ± 0.3 GeV, for RHIC and LHC data, respec- 
tively. The values of n are, 0.26 ± 0.01 and 0.01 ± 0.00, 
for RHIC and LHC data, respectively. Figure 0] (b) 
shows the 0.5 A^pa^-normalized ^^ for Pb+Pb col- 
lisions at ^/snn = 2.76 TeV for LHC and for Au+Au 
collisions at ^/snn = 200 GeV for RHIC as a func- 
tion of collision centrality. A similar function as is 
given by Eqn. \5\ has been fitted to both the data sets. 
The fitting parameter n, has the values —0.08 ± 0.01 
and -0.07 ± 0.01 for RHIC and LHC data sets, respec- 
tively. The parameter B, has the values 5.96+0.34 GeV 
and 6.01 ± 0.14 GeV for RHIC and LHC data, respec- 
tively. With top peripheral collisions as exceptions, the 
0.5 A^ q _p ar t-normalized ^ both for RHIC and LHC 
data within errors, show a centrality independent num- 
ber of constituent quarks scaling behavior. The value 
of the power, n in the power law gets smaller while go- 
ing from A^N-part to Af q _ par t-normalization, thereby mak- 
ing the spectra in Af q _ pa rt natter compared to A^-part- 
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FIG. 5: (Color online) AN-part and A q _ par t normalized —^ L 
as a function of collision energy from lower AGS, SPS and 
RHIC to LHC energy. The collision data are compared with 
IP-saturation model and hybrid function estimations. 



normalization. The shapes of the centrality dependent 
A^N-part or Ar q _ part -normalized ^- and ^ show a 
universal feature. The relative enhancement of iVj^— part- 
normalized ^^ in A+A collisions w.r.t. that in p + p 
collisions is the same at RHIC and LHC energies. This 
has been explicitly seen in the framework of models based 
on multiple scattering and mutual boosting of saturation 
scales and shadowing (l9| |. 

In Figure ^j^ normalized to both A^-part and 
■Aq—part is shown as a function of collision energy. Up 
to top RHIC energy Af N _ part -normalized —jf*- is well de- 
scribed by logarithmic function given by Eq. [5] How- 
ever, this function fails to describe the ALICE data for 
2.76 TeV. The estimation based on logarithmic func- 
tion for the Af N _part-normalized ^^ is 1180 ± 22 for 
y/sNN = 2.76 TeV, which underestimates the data by 
26%. A power law function (dash-dotted line) given by 
Eq. discribes both RHIC and LHC heavy-ion data 
but overestimates lower-energy measurements. Contrary, 
the collision energy dependence of £ ^ k normalized by 



q_part is very well described by logarithmic function for 
the whole range of energies under discussion with little 
deviation towards LHC energies. In this case, the power- 
law gives a better description of data at the level of quark 
participants. The power in the power law function de- 
creases from A^N.part to iV q _ par t normalization, thereby 
going towards a flatter behavior as a function of collision 
energy. The predicted value of —^ for Pb+Pb collisions 
at ^/snn = 5.5 TeV based on the extrapolation of power- 
law function fitted to A/q_ par t and A^-part-normalization 
is around 2116 ± 52 and 2005 ± 28, respectively. On the 
other hand, the value predicted based on a logarithmic 
extrapolation of AN-part-normalization data is 1310+22. 



Looking at the low-energy and high-energy behaviors of 
charge particle production being well described by a log- 
arithmic function and power-law functions, respectively, 
we have tried to fit a hybrid function (a combination 
of both) given by Eq. [8] Both A^-part to iV q _part- 
normalized , ch for the whole range of energies is very 
well described by this hybrid function, as is seen from 
Fig. [SJ At the level of nucleon and quark participants 
this hybrid type of function describes the charge parti- 
cle production universally. At the level of Af q _ par t and 
A^N-part-normalization, the ^j^ estimated by the hybrid 
function for Pb+Pb collisions at ^/snn = 2.76 TeV de- 
viates from the experimental data by 0.06% and 0.44%, 
showing an excellent agreement with the data. The pre- 
dicted values of ^^ for Pb+Pb collisions at ,/snn = 5.5 



drj 

TeV based on the extrapolation of the hybrid function 
fitted to A^-part and A^-part-normalizations are around 
2836 ± 85 and 2168 ± 73, respectively. The predictions 
from IP-saturation model for the top RHIC energy and 
higher, are also shown for a comparison with the corre- 
sponding experimental data. 
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FIG. 6: (Color online) Predictions for mid-rapidity d ch for 
central Pb+Pb collisions at -^/snn = 5.5 TeV at the LHC. On 
the left the names of the authors and on the right, the observ- 
ables and the centrality definitions are shown. The error bars 
show the uncertainties in the predictions. Our predictions 
both from nucleon and quark participant normalizations are 
shown as filled triangles. 



y = C(s NN ) n (7) 



y = A + B In ^/IWn + C(y/s^) n (8) 

The form of the above hybrid function is motivated by 
the energy dependent behavior of charged particle pro- 
duction and in addition, by the following results from 
the non-equilibrium statistical relativistic diffusion model 
(RDM) [2fJ. It is explicitly shown by the works of G. 
Wolschin et al. [H| that the RHIC and LHC multiplicity 
could be explained by the combination of a mid-rapidity 
gluonic source and a fragmentation source. Mid-rapidity 
gluonic sources predict a power-law type behavior of 
charged particle multiplicity, whereas the fragmentation 
sources predict a logarithmic behavior. Energies lower to 
that of RHIC, where particle multiplicity seems to obey 
a logarithmic behavior, could easily be explained by frag- 
mentation sources. In Figure [SJ the predictions for mid- 
rapidity -j^ for central Pb+Pb collisions at ^/snn = 5.5 
TeV at the LHC, are shown and compared with our pre- 
dictions from the hybrid function at nucleon and quark 
participant level. Other data are taken from Ref. [22j. 

A similar analysis is done for =§Z- -normalized to 
A^-part an d Af q _ part and is plotted as a function of colli- 
sion energy in Figurc[7l Up to top RHIC energy A^-part- 
normalized ^f^ is well described by logarithmic function 
given by Eq. [6] However, this fails to describe the AL- 
ICE data for 2.76 TeV. The estimation of ^ based on 



logarithmic function for the A^N-part-normalized ^p is 
1114 + 24 GeV, which underestimates the data by 47% 
for ^/snn = 2.76 TeV. A power law function (dash-dotted 



line) given by Eq. [7J discribes both RHIC and LHC 
heavy-ion data but underestimates lower-energy mea- 
surements. Contrary, the collision energy dependence 
of ^j^ normalized by A7q_ par t is very well described by 
logarithmic function up to top RHIC energy and shows 
considerable deviation at energies higher to top RHIC 
energy. The power law function given by Eq. [Jj de- 
scribes the A'q-part-normalized ^j 21 data quite well. Like 

^j^ 1 - , the power in the power law function, decreases from 
A^-part to iVq-part normalization, thereby going towards 
a flatter behavior as a function of collision energy. The 
most interesting observation is that the power-law mo- 
tivated function better describes the transverse energy 
production with collision energy compared to the loga- 
rithmic function. The predicted value of ^^- for Pb+Pb 
collisions at ^/snn = 5.5 TeV based on the extrapola- 
tion of power-law function fitted to Af q _ par t and iVji— part- 
normalized #£ arc around 2683 + 62 GeV and 2760 + 50 

ari 

GeV, respectively. Like the case of charge particle pro- 
duction, looking at the low-energy and high-energy be- 
haviors of the transverse energy production being well 
described by a logarithmic function and power-law func- 
tions, respectively, we have tried to fit the same hy- 
brid function given by Eq. [S] The A^-part-normalized 
^j 2 - for the whole range of energies is very well de- 
scribed by this function, as is seen from Figure [7] The 
A^q-part-normalizcd ^^ is also very well described by 
the hybrid function. Irrespective of nucleons or quarks 
are the sources of particle production, this hybrid func- 
tion describes the charged particle and transverse energy 
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FIG. 7: (Color online) iVN-part and -/Vq_ par t normalized -^ s - 
as a function of collision energy from lower AGS, SPS and 
RHIC to LHC energy. The collision data are compared with 
the EKRT model and the hybrid function predictions. 



production for all range of energies starting from lower 
AGS, SPS, RHIC to LHC. The estimations of ^§f- from 



part 

2.76 



the hybrid function fitting to the -/V q _ part and JVn 
normalized data for Pb+Pb collisions at ^/snn = 
TeV deviate from the CMS collaboration experimental 
data (H by 0.05% and 1.18%, respectively The pre- 
dicted value of ^j^ for Pb+Pb collisions at ^/snn = 5.5 
TeV based on the extrapolation of the hybrid function 



lisions. This quantity depends on the initial state of the 
collision and the viscosity of the matter as it expands to 
its final state, when it is observed by the detectors. This 
observable when studied as a function of collision energy, 
shows two regions of interest. The first one from the 
lower SIS energies to SPS energies shows a steep increase 
of Et/NcIx values, thereby indicating that the mean en- 
ergy of the system increases. In the second region, from 
SPS to top RHIC energy, £+/V c h shows a very weak col- 
lision energy dependence, i.e. like a saturation behavior. 
In this region the mean energy doesn't increase, whereas 
the collision energy increases. This indicates that the in- 
crease in collision energy results in new particle produc- 
tion in this domain of energy. This behavior has been 
well described in the context of a statistical hadron gas 
model (SHGM) [H. In the framework of SHGM, it has 
been predicted that E^/N c ^ would saturate at energies 
higher to that of top RHIC energy with a limiting value 
of 0.83 GeV [2y|. This goes inline with Hagedorn's con- 
jecture of a limiting temperature (Th = 165 MeV) for 
elementary and nuclear collisions at higher collision en- 
ergies [28| and the subsequent observation of a saturation 
behavior in the freeze-out temperature (Tf ) in heavy- ion 
collisions 0. The SHGM prediction of E T /N ch at LHC 
collision energy of 2.76 TeV has been found to be incon- 
sistent with the experimental observation of its value, 
which is 1.25 ± 0.08 GeV [23J|. This could be because of 
more high-px particle and jet production at LHC, which 
makes LHC measurements quite different from that of 
RHIC. The estimated values of Et/N^ from the hybrid 
equation fitting to the corresponding Vq-part-normalized 
data, for Pb+Pb collisions at 2.76 and 5.5 TeV collisions 
energies are 1.33 ± 0.05 and 1.49 ± 0.05 GeV, respec- 



fittcd to A r a _ Dart -normalization is around 3056 + 44 GeV. ,. , n , , -, ,, AT ,. , , , 

_ q part , . tively. Here we have taken the A* a _ Dar t-normalizcd data 



The transverse energy production as a function of en 
ergy in the framework of EKRT model, which is based 
on initial state gluon saturation [l8[ , is given by, 



dE T 
dr\ 



0.46 A°- 92 {^) 0A0 [l - 0.012 In A + 0.061 lny/a\, 

(9) 

where A is the mass number of the colliding nuclei for 
a symmetric collision species. This takes into account a 
reduced initial number of scattering centers in the nuclear 
parton distribution functions. Here we have compared 
the EKRT model estimations for transverse energy for 
all collision energies under discussion. It is observed that 
the transverse energy production at higher energies is 
very well described by EKRT model. However, EKRT 
model fails to describe the low energy data. This is an 
indication of gluon saturation effect at higher collision 
energies. 

The ratio of pscudorapidity densities of transverse 
energy and number of charged particles i.e. Et/Nc^ 
has been studied both experimentally [23|-[25j and phe- 
nomenologically [2^, [27| to understand the underly- 
ing particle production mechanism. This observable is 
known as global barometric measure of the internal pres- 
sure in the ultra-dense matter produced in heavy-ion col- 



and the corresponding fitting of the hybrid function, as 
it describes both charged particle and transverse energy. 
This shows that the E^/N^ shows a very strong depen- 
dence on collision energy beyond the top RHIC energy, 
thereby creating a third region in the Et/N^ Vs colli- 
sion energy diagram. If the freeze-out occurs on a fixed 
isotherm for all energies, it is expected that the energy 
per particle will be the same. A comparison of Tf and 
Et/N^ versus collision energy indicates that from SPS 
to top RHIC energies the T{ and the value of Et/N^ 
are almost independent of collision energy. Beyond top 
RHIC energy, the value of Et/N^ gets a jump, which 
indicates that the T{ at LHC should be different than 
that is at RHIC. We have estimated the Tf for the LHC 
energies from the SHGM parametrized equations used in 
Ref. j29[, which are given by 



T(u B ) = a-bii 2 B - c/j.%, 



(10) 



where a = 0.166 + 0.002 GeV, b = 0.139 + 0.016 GeV~ 
and c = 0.053 ± 0.021 GeV 3 and 



I^b(Vs) 



l + e^/Ts 



where d 



1.308 ± 0.028 GeV and e 



(11) 



0.273 ± 



0.008 GeV 1 . The value of the freeze-out temperature, 
T fo , for RHIC 200 GeV, LHC 2.76 TeV and 5.5 TeV are 
around 166 ±2 McV. However, the T{ obtained from the 
SHGM estimation of fitting to experimentally obtained 
particle ratios at RHIC and LHC is T io = 164 MeV, with 
fi B is 24 MeV and 1 MeV, respectively [H ED]. The 
SHGM predictions for Tf is consistent with its measure- 
ments at LHC. Contrary, the SHGM predicted value of 
Et/Nc^, which is 0.83 GeV is inconsistent with the cor- 
responding measurements at LHC, thereby indicating a 
new physics governing particle production and thermal- 
ization at LHC. The observation that T{ doesn't change 
while going from RHIC to LHC energies, when Et/N^ 
gets a jump, needs to be understood from thermodynam- 
ics view point. Tf being independent of collision species 
and energy is a consequence of the fact that the system 
produced in heavy-ion collisions irrespective of the initial 
conditions comes to the same final state. Et/N^, the 
mean energy per particle, getting enhanced beyond top 
RHIC energy, could be a consequence of LHC being com- 
pletely dominated by high-px jets and as a barometric pa- 
rameter, its increase implies, higher initial pressure being 
created in LHC collisions compared to that of RHIC. This 
also is reflected from the values of radial flow parameters 
measured at RHIC (< &v >= 0.59 ± 0.05 [H) and LHC 
(< j3 T >= 0.65 ± 0.02 [3l|). Radial flow for central colli- 
sions at LHC is 10% higher than the corresponding value 
at RHIC in the same j?t range. In addition, this rise in 
Et/N^ could be understood from the gluon saturation 
picture as follows. The number of gluons inside the col- 
liding particles increases with increase in collision energy. 
At very high energy, the gluon creation and annihilation 
balance out leading to a saturation in gluon number. As 
gluon density does not have a linear collision energy de- 
pendence, the saturation scale depends non-linearly on 
the center of mass energy i.e. Et/N^ = fc(y / s) A , where 
k and A are constants [33j. This function when fitted to 
the top RHIC and LHC data, where a gluon saturation 
is expected, describes the energy dependent behavior of 
Et/N,^ quite well. The value of A obtained from this 
fitting is found to be 0.12, which lies between A = 0.15 
(IP-sat mo del p rediction) and A = 0.11 (b-CGC model 
prediction) [3j,[3J]. Hence, gluon saturation picture at 
LHC energies may provide an explanation to the jump 
in E T /N ch [H. 



IV. SUMMARY AND CONCLUSION 

In this paper, we have estimated the number of con- 
stituent quarks participating in a nucleus-nucleus colli- 
sion in the framework of a nuclear overlap model. The 
overlap model estimation of number of nucleon partici- 
pants agrees well with the corresponding estimation from 
ALICE experiment for Pb+Pb collisions at -y/SNN = 2.76 
TeV. The ratio of number of quark participants and num- 
ber of nucleon participants shows a non-linear increase 
with collision centrality for a particular collision energy. 
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FIG. 8: (Color online) The ratio of —^- and d ch as a func- 
tion of center of mass energy. Experimental data are com- 
pared to the predictions from thermal model, gluon satura- 
tion model and the estimations obtained in the framework 
of the hybrid model fitting to transverse energy and charged 
particle data. 



The centrality dependence of iVN- pa rt-normalizcd ^f^ 
both for Au+Au collisions at -^/snn = 200 GeV and 
Pb+Pb collisions at ^/snn = 2.76 TeV, which show a 
increase with centrality, turns out to be a centrality in- 
dependent scaling behavior when normalized to A q _ part . 
This gives the indication of partonic activities at RHIC 
and LHC energies. Similar is the observation for =§£• 
both for the RHIC and LHC data. The charged parti- 
cle production as a function of collision energy has been 
studied both in the nucleon and quark participant frame- 
work. A logarithmic function describes the low-energy 
production of charged particles quite well up to top RHIC 
energy but fails to describe the LHC data in the nucleon 
participant sector. However, in the quark participant 
sector, the logarithmic function describes the data well 
within errors for energies up to LHC collision energy of 
2.76 TeV. A power-law function however, describes the 
high energy data with an excellent agreement with LHC 
2.76 TeV measurement of — f^. On the other hand, the 
power-law function underestimates the low energy data. 
A hybrid function, which is a combination of logarithmic 
and power-law in collision energy, describes the particle 
production at all energies extremely well both for the nu- 
cleon and quark participant sector, indicating a universal 
law of particle production. We have given the predictions 



for 



dN cl 
di] 



for Pb+Pb collisions at 



5.5 TeV from 



this universal hybrid function. The hybrid function de- 
scribes the transverse energy production at all energies 
under discussion both in the nucleon and quark sectors. 
However, the EKRT model based on a initial state gluon 
saturation and the power-law function better describes 
the high energy behavior of transverse energy produc- 



tion. This indicates the gluon saturation effects at en- 
ergies higher to RHIC. We have estimated the vaiue of 
^t for Pb+Pb collisions at y/s^ = 2.76 TeV by us- 
ing this hybrid model which agrees with the correspond- 
ing LHC measurements. And thus, we give a predic- 
tion of ^j^ for the Pb+Pb collisions at ^/snn = 5.5 
TeV. The estimated value of E T /N ch for Pb+Pb col- 
lisions at -y/sNN = 2.76 TeV estimated by the hybrid 
function agrees with the corresponding measurement by 
CMS collaboration at LHC, showing a jump in its en- 
ergy dependence behavior. This indicates a higher initial 
pressure of the ultra-dense matter at LHC compared to 
that is formed in RHIC energies. In this framework, we 
give the prediction for this barometric observable for the 
future LHC Pb+Pb collisions at ^/sii = 5.5 TeV. We 
have seen explicitly that the discussed hybrid function 
also describes the strange and multi-strange particle pro- 
duction at all energies. In addition, we have also seen 
explicitly that this function describes the particle pro- 



duction in elementary e + e~ , p + p(p), e p, p + p and 
p + A collisions quite well for all energies. In conclu- 
sion, the charged particle particle and the transverse en- 
ergy production mechanisms are very well described by 
the hybrid function for all range of energies and collision 
species starting from elementary to heavy-ion collisions, 
indicating a universal law of particle production mecha- 
nism. However, the observed value of E^/N^ at LHC, 
which is not explained by SHGM, is well reproduced in 
a gluon saturation picture. 
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